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Abstract — In this paper, we investigate the pilot position selec- 
tion (PPS) problem in the frequency domain pilot multiplexing 
technique (FDPMT) for channel estimation of single-carrier block 
transmission with frequency domain equalization (SC-FDE). Un- 
like the widely accepted point of view that the conventional PPS 
technique is the optimal one, this paper for the first time questions 
this viewpoint and demonstrates the suboptimal property of the 
conventional PPS technique. To further improve the performance 
of the FDPMT, an optimal PPS technique is proposed based on 
the minimization of the average bit error rate (BER). Compared 
with the conventioanl PPS technique, the proposed PPS technique 
exhibits better performance with similar complexity. 

I. Introduction 

Unlike orthogonal frequency division multiplexing (OFDM) 
systems, the SC-FDE transmits the data block-wisely in the 
time domain while performs the channel estimation and equal- 
ization in the frequency domain. Compared with the OFDM, 
the SC-FDE has similar overall implementation complexity 
but exhibits superior system performance [1]. To improve 
the data transmission performance, coherent demodulation is 
commonly applied for both OFDM and SC-FDE systems. The 
successful application of such technique critically depends on 
the accuracy of the channel estimation. For the past years, 
plenty of investigations have preferred utilizing the frequency 
domain multiplexed (FDM) pilots insertion technique for the 
channel estimation of OFDM systems while preferred using 
the time domain multiplexed (TDM) pilots insertion tech- 
niques for the channel estimation of SC-FDE systems [1][2]. 

In [3], two types of FDM pilot insertion techniques, i.e., 
frequency expanding technique (FET) and frequency domain 
superimposed pilot technique (FDSPT), were proposed for 
the first time for SC-FDE systems to track the channel state 
information (CSI). Unlike the FET, the FDSPT does not 
need extra spectrum for the pilot tones. Instead, the equally 
spaced pilots are superimposed on the data tones starting from 
the first tone for all SC-FDE blocks and thus the FDSPT 
keeps the spectral efficiency as the original SC-FDE systems. 
However, the simulation results in [3] [4] have shown that the 
introduction of the mutual interference between data symbols 
and pilots in the FDSPT degrades the BER performance. 

To improve the performance of FDSPT, a new FDPMT 
was proposed in [5]. Unlike the FDSPT, the FDPMT first 



eliminates some equally spaced data tones whose initial po- 
sition is selected based on the proposed PPS technique at 
the transmitter and then inserts pilots block by block. Based 
on such PPS technique and by iteratively reconstructing the 
distorted data symbols at the receiver, the FDPMT provides 
better BER performance than the FDSPT and even approaches 
the lower bound of SC-FDE systems after many iterations. 
Note that compared with the FDSPT, the essential innovation 
of the FDPMT is the development of the PPS technique. Up 
to the present, it has been widely accepted in academia [6]- 
[10] that the proposed PPS technique in [5] is the optimal one 
in achieving the ultimate system performance. 

In this paper, we reinvestigate the principle of the FDPMT 
and for the first time question the aforementioned widely 
accepted viewpoint based on our observation. Specifically, 
different from the superimposed pilot scheme in an OFDM 
system where the distortion caused by the insertion of pilots 
can be treated as noise [11], the loss of data tones in the 
superimposed pilot scheme for a SC-FDE system actually 
results in the deterministic interference. This is because that 
the decision for SC-FDE signals is performed in the time 
domain. The above mentioned observation indicates that the 
current PPS technique based on the minimization of the 
distortion of SC-FDE signals [5] would not be equivalent 
to the optimization for the symbol decision. To prove the 
suboptimal property of the PPS technique in [5], 16-QAM 
modulation is adopted as an example for the derivation and 
analysis throughout this paper. 

The suboptimal property of the current PPS technique 
motivates us to propose an optimal one. Therefore, in this 
paper, an optimal PPS technique is proposed, where the 
selection of the initial position of pilots is based on the 
minimization of the average BER in each SC-FDE block with 
pilots. Compared with the significant performance gain of the 
proposed PPS technique over the current one, the increase of 
the complexity of the proposed PPS technique over the current 
one is marginal. Moreover, the proposed PPS technique with 
no iteration at the receiver even expresses better performance 
than the current one with one iteration for higher signal-to- 
noise (SNR). This means that for higher SNR the proposed 
PPS technique can make the SC-FDE system express similar 
complexity to the system applying the current PPS technique 
while still preserve the performance gain. 
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II. System Models 

Fig. 1 depicts the block diagram of the SC-FDE transmitter 
with FDM pilots [5][8]. The SC-FDE transmit block is of 
length TV and denoted by s = [sq, . . . , sjy-i] with each symbol 
drawn from a complex valued alphabet. The TV-point discrete 
Fourier transform (DFT) of s is given by 
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For the FDPMT, some data tones are deliberately eliminated 
to insert pilots in substitution, and thus the actual frequency 
domain transmit sequence becomes 
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where C = [Co, . . . , Cn p -i\ is the pilot sequence of length 
TV p , M = N/N p is the spacing between two adjacent pi- 
lot tones, and do is the initial pilot position. By defining 
X m = [So, • • • j S m -i, 0, . . . , SM+m-ii Oj • • • , <5iv-i] with 
m € [0, M) denoting the initial pilot position and its inverse 
DFT (IDFT) as x' m = [x' mi o, ■ ■ ■ ,x' mj N-i], the distortion 
level due to the removal of the data tones at indexes t/j m = 
{m, M + m, . . . , (TV p — 1) M + m} is quantified in [5][8] as 



= E 



N-l 



, and accordingly d a is selected 



based on its minimization as 



da = argmine m = argmin$^$, 
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where # m = [S m , S M+m , . . . , S, 



] and (•) de 



(N p -l)M+m 

notes complex conjugate transpose. It is easy to find that the 
above strategy is in fact equivalent to minimize the loss of the 
power of s since we have Y^nZo I s " 1 2 = ^k=o \Sk^ based 
on the Parseval theorem. Before transmission, Xd 0 is first fed 
into an TV-point IDFT and then appended by the length-G 
cyclic prefix (CP). 

At the receiver, an TV-point DFT is performed to convert 
the received signal without the CP to the frequency domain as 
Rrf 0 = [Rdofl, ■ ■ ■ ,Rd 0 ,N-i}- By adopting the pilot position 
detection techniques proposed in [8]-[10], the pilot tones, 
i.e., {Rd 0 ,k, k (mod M) = do}, can be easily and properly 
extracted from R<j 0 . Following that, the channel frequency 
responses (CFRs) on the pilot tones are estimated by dividing 
the extracted signals by the prior known pilots and those on 
the data tones are then obtained by performing interpolation. 
Many interpolation techniques can be used, such as linear 
interpolation, trigonometric interpolation, and etc., and their 
performance comparison has been investigated in [12]. For 
analytical simplicity, in this paper we assume perfect CSI. 

Once CSI is known, the decision of the SC-FDE symbols 
should be made coupled with the iterative signal reconstruction 
(ISR) technique to improve the overall system performance [8] 
[13]. This is implemented as follows. Firstly, by equalizing the 
elements on data positions and nulling those on pilot positions, 
the outputs in the frequency domain can be readily obtained 
as 
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Fig. 1 . The block diagram of the transmitter with the FDPMT. 



with k = 0, 



E 



TV — 1, where Hk is the CFR on the fcth tone, 
= 1, and Vk is the zero-mean complex additive 

white Gaussian noise (AWGN) with variance a^. Without loss 
of generality, zero forcing (ZF) equalization is assumed in (4). 
Secondly, the initial estimates of SC-FDE symbols are derived 
by performing coherent demodulation to the TV-point IDFT of 
those outputs in (4). Finally, the ith (i > 1) estimates are 
iteratively obtained by the IDFT of 
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with k = 0, ...,TV - 1, where {S^'} 1 ^ are the DFT 

of the (i — l)th estimates, i.e., {s^rTn ln=o- Notable is that 
by jointly considering the complexity and performance, the 
suggested iteration number should not exceed one. 

III. OPTIMAL PILOT POSITION ANALYSIS 

The non-linear operation of ISR would easily incur error 
propagation, consequently to guarantee the system perfor- 
mance, reliable initial estimates must be satisfied [13]. To 
get further insight, let us consider the initial input for a 
demodulation, which turns out to be the IDFT of (4), i.e., 
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0 for k (modM) 



distortion caused by superimposing pilots, and v'd 0 , n = 
T^F Ysk^o 1 ( a d 0 ,kVk/H k ) e^ 2 ^ denotes the equalized noise. 

Intuitively, the PPS criterion according to the minimization 
of the power loss in (3) should be significantly helpful in 
reconstructing the desired signals {s^nTJo 1 wnen {s'do^l^Jo 1 
in (6) are considered as part of noise. However, for SC-FDE 
systems in which {shI^Tq 1 are drawn from a finite alphabet 
set, the distortion actually contributes as a deterministic inter- 
ference from coherent demodulation perspective. This raises 
several questions, such as whether the minimization of the 
distortion of SC-FDE signals optimizes the symbol decision, 
and if not how to design the PPS technique to obtain the 
optimal symbol decision. In this section, we will address the 
above questions. 
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Fig. 2. Signal space diagram for 4-PAM. 



A. Description of the suboptimal property of the current PPS 
technique 

Without loss of generality, suppose that 16-QAM modula- 
tion is adopted. In this case, a modulated symbol consists of 
4 bits, i.e., i\iiq\qi, with the first two bits representing the 
in-phase (I-) component and the last two bits representing the 
quadrature (Q-) component. Since the analysis for I-channel 
and Q-channel scenarios is similar, the I-channel scenario is 
taken as an example for detailed analysis. Fig. 2 shows the I- 
channel signal constellation (4-PAM) using Gray mapping and 
its decision boundaries, where l\ and 12 denote the regions in 
which i\ = 1 and i% = 1, respectively, and d = with 
Eb denoting the signal energy per information bit. 

To verify that the PPS criterion in (3) is a suboptimal 
strategy, we resort to the following example. 

Example: Suppose that two options are available for su- 
perimposing pilots: one is at tones ip mi and the other is at 
tones il' m2 - To make this example easy to follow, we assume 
the accumulated distortion level on all SC-FDE symbols 
excluding the real part of the nth sample in one block for 
both options are equal, i.e., Yl u ?n l s 'mi,«| 2 + |3 [s'™,,,,]! 2 = 
J2 u ^n \ s 'm 2 ,u\ 2 + |3 [s'm 2 ,n]| 2 - In this case, the selection 
between these two options mainly depends on the distortion 
level on the real part of the nth sample, s r n . If s r n = 3d and 
s mi,n < 0 < s' r m2 n (s' r X n = U[s' x ,n\ , X G {mi,m 2 }), 
i.e., the signal point labeled "10" of the nth sample is shifted 
to the left for the first option and to the right for the second 
option due to the distortion. It can be concluded that compared 
with the unshifted case, the error probabilities of both bits, 
i\ and 12, would become bigger for the first option while 
smaller for the second option. This is because the signal 
point gets closer to the decision boundaries or even moves 
to the wrong decision range for the left-shifted case while 
gets further from the decision boundaries for the right-shifted 
case. However, when the left-shifted amount is smaller than 



the right-shifted amount, i.e., s 



< s' 



7712 ,71 I 



the first 



option will be chosen instead according to the current PPS 
criterion. Similarly, for s r n = —3d and s' r min < 0 < s' r m2 „, 
it will be seen that the current PPS criterion would erroneously 
pick up the second option provided that the right-shifted 
amount is smaller than the left-shifted amount. The analysis 
for the suboptimal property of the current PPS criterion in the 
case of s r n = ±d and s' r min < 0 < s' r m2n seems a bit tedious 
since the comparison of the error probabilities for both options 
is highly relevant to the specific realization of channel, we omit 



X 771,71 ~T" ^ 771,77 5 Where V 777,77 

~ —o^ due to a linear 



it here due to the space limit. Note that there exist many other 
counterexamples for the validation of the suboptimal property 
of the current PPS technique. 

B. An optimal PPS technique 

Subsection III. A clearly demonstrates that the minimization 
of the distortion of SC-FDE signals cannot optimize the 
symbol decision. This subsection will give the answer about 
how to design the PPS technique to obtain the optimal symbol 
decision. From the above analysis, it is clear that an optimal 
PPS criterion should aim at minimizing the overall impact of 
the superimposition pilots on the SC-FDE systems, which can 
be characterized by using the average BER in one SC-FDE 
block. Motivated by this, in the following we first derive the 
closed-form expression of the average BER and then propose 
the realization of the optimal PPS criterion. 

Let us begin with the simple case of an AWGN channel. 
Without loss of generality, consider the error probability 
of the nth sample, s n , with the assumption that pilots are 
inserted at tones ip m . The initial input for demodulation can 
be obtained from (6) as §m]n 
is AWGN with variance a%, 
combination of Gaussian variables. For the I-component of 
s n, s r n = Sft [s„], there are two possible cases (Case I and 
Case II) for the probabilities, P^ n , that the bit i t (I = 1,2) 
is in error. To simplify the following analysis, we introduce 
parameters S, n ,ia (I = 1,2) to further fractionize each 

class into four situations according to the possible positions 
of and x' r mn {x' r mn = 3? [a/ m ,„]) in the constellation with 
each corresponding to {;„,/,! = £„,;, 2 = 1, fn.i.i = -1 
and £ n> i t 2 = 1, £ n ,l,i = €n,l,2 = -1, and £„,; ; i = 1 and 
£n,',2 = — 1, which are in sequence referred to as "si", "s2", 
"s3", and "s4", respectively, in this paper. Here, £ n ,i,i = 1 
indicates that s r n belongs to any signal point in the region %i 
otherwise £ n ,i,i = — 1 and £ n ,i,2 = 1 stands for x' r m n lying 
in the region ij_ otherwise £ n ,',2 = — 1- 

For Case I, only the bit i\ is considered. The constellation, 
as shown in Fig. 2, can be separated into two regions based on 
the decision boundary represented by the solid line at the ori- 
gin. As mentioned above, in total four situations for Cass I will 
be encountered. Specifically, for "si", i.e., £ n ,i,i = £,1,1,2 = 1, 
a bit error will occur if the noise is smaller than \—x' T m1 \. 
In this case, the probability that the bit i\ for "si" is in error 
is [14] 



JV 



^777,71 ^1 




(7) 



where Q (x) = (1/V2tt) f™ e^^dt. For "s2", i.e., £„, M = 
— 1 and £71,1,2 = 1, a bit error will occur if the noise exceeds 
(~ x m,n). yielding 



p7l I 
1 777,7l|s2 



= l-Q 



2 (^'777,71) 



(8) 



Due to symmetry, the bit error probability of i\ for "s3", 
i-e., 67,1,1 = £n,i,2 = -1. and "s4", i.e., £„,i,i = 1 



5068 



IEEE ICC 2014 - Wireless Communications Symposium 



and £„.i,2 = — 1, will be equal to that for "si" and "s2" 
respectively, i.e., P^J s3 = P£,J a i and P^„U = P&, n \«2. 
Therefore, putting these four situations together, we can obtain 
a generic form for the error probability of the bit i\ as 



pii 



1 — 



6n,l,l.&i 




Note that there are much more possible manners to generalize 
the form. For example, we could have used £n,i,i/£n,i,2 

or (|£n,i,i + £n,i,a| ~ !) in P lace of Cn,i,i^n,i,2, where | ■ | 
denotes the absolute value. 

For Case II, we consider only the bit i% while ignore the 
bit i\. The dashed lines crossing —2d and 2d are given as the 
decision boundaries of the decision region in Fig. 2. Firstly, 
let us focus on "si", i.e., £„,2,i = £ n ,2.2 = 1. A bit error will 
occur if the noise is smaller than (—2d — x' r m n ) or exceeds 
(2d — x' r m n ). Therefore, the error probability of the bit i 2 for 
"si" can be expressed as 



Pm,n\sl — Q 



l2(x' r mtn -2d)_\ Q (j2(x' r mtn + 2dY 



(10) 



Secondly for "s2", i.e., £ n ,2,i = —1 and £n,2,2 = 1, a bit 
error will occur if the noise is larger than {—2d — x' r m n ) and 
smaller than {2d — x' r m „), resulting in 



PZn\s2 = l~Q 




Thirdly in a similar manner, we can calculate the bit error 
probabilities of i 2 for "s3", i.e., £ n ,2,i = £n,2,2 = — 1, and 
"s4", i.e., i n 21 = 1 and £ n 2.2 = -1, as 



P 12 I i 

1 m,n\si 



and 



P 22 I A — 1 

r m,n|s4 — 1 




'2(^„-2d) z \ _ Q //2( a /^ >n +2d) i 



(13) 

respectively. Finally, careful inspection of (10)-(13) reveals a 
generic expression for the error probability of the bit i 2 , which 
is of a form 



Since the demodulation of two quadrature channels are 
independent, from symmetry the error probabilities of the bit, 
qi, and the bit, q 2 , P^ n and P$ n , can be calculated similar 
to (9) and (14) but with the terms x' r m n , £ n ,M> and £ n ,i,2 
be substituted by x' l m n , ( n ,is, and Cn,i,2, respectively, where 
C«,i,i and (n,L2 have similar definitions to £ n ,j,i and £„,; |2 , 
respectively (1 = 1, 2), and x n mn = 3 [x'm.n]- It can be seen 
from (9) and (14) that the expression of bit error probability 
consists of a sum of Gaussian g-Function, with a general form 
Q (y/2p), where p depends on SNR, distortion, and decision 
boundary. 

In the case of a flat fading channel, the CFRs for all tones 
are the same, i.e., Hk = H for all k = 0, . . . , N — 1. The error 
probabilities of bits i\ and i 2 carried on s r n can be obtained 
by first substituting the common term in (9) and (14) with 
<j%i/\H\ 2 and then taking the expectation with respect to \H\ 2 . 
For example, assume a Rayleigh fading channel in which \H\ 2 
is characterized by exponential distribution. Utilizing the result 
J 0 °° Q (v / 2m/5) e~ u du = | ^1 — 'jf^, the expression for 
the error probability of each bit for a Rayleigh flat fading 
channel can be derived by using the form | ^1 — J m 

place of the form Q (\/2p) in an AWGN case. The case for 
bits qi and q 2 is straightforward. 

In the case of a frequency-selective fading channel, the 
conditional output SNR of Sm]n on the CFRs can be readily 
expressed from (6) by 



/JV-l 2 

7 =4E h ( V am ^ v ' 

Vfe=0 1^*1 



(15) 




which is the harmonic mean of the SNRs on data tones. In- 
terestingly, this SNR expression has been also encountered in 
the multi-hop CSI-assisted amplify-and-forward (AF) relaying 
[15][16]. As demonstrated in [17], it is difficult to derive 
the analytical BER by adopting the commonly used moment- 
generating function (MGF) method in [14] since there is so far 
not closed-form expression for the MGF of the instantaneous 
SNR as in (15) except for two chi-squared random variables, 
i.e., N = 2 [18]. The authors in [17]-[20] have suggested an 
alternative but mathematically tractable method based on the 
characterization of the statistical properties of the equalized 
noise, i.e., v' m , n . However, this method is only applicable 
under the assumption of independence among the CFRs and 
not straightforward for dependent case. To cope with this 
problem, we notice that it is inessential to obtain the explicit 
BER while feasible to use its tight lower bound instead for 
the comparison of the achievable performance for different 
initial pilot positions. By adopting the tight upper bound on 
the output SNR proposed in [21] that is the minimal value 
among the SNRs on the data tones (and thus irrelevant to 
the correlation of the CFRs), the lower bounded BER (or say 
approximated BER) for a frequency-selective fading channel 
can be readily obtained using the results given in fiat fading 
channel case. The feasibility of this treatment will be verified 
by the simulation in Section IV. 

Up to the present, the error probability of each bit carried 
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Fig. 3. BER performance comparison between the current PPS and proposed 
PPS techniques in an AWGN channel. 




SNR (dB) 

Fig. 4. BER performance comparison between the current PPS and proposed 
PPS techniques in a Rayleigh flat-fading channel. 



on s n has been derived. Based on this, the average BER in 
one block due to the removal of the data tones at index ip m 
can be written as 

1 N ^ 

n=0 

By following steps similar to those given above, the average 
BER can be easily obtained for higher-order QAM, e.g., 64- 
QAM, 128-QAM, and so on. 

According to current analysis, the optimal initial pilot 
position should be selected to minimize the average BER, such 
that we have 

d 0 = argminP mje . (17) 

m 

IV. Simulation Results and Discussions 

In this section, BER simulations are provided to verify the 
utility of the proposed PPS technique. An uncoded SC-FDE 
system with 16-QAM modulation and FDPMT is considered. 
The transmit block consists of 512 modulated symbols and 
a guard of length 26. The pilot overhead ratio is chosen 
as 3.125% (N p = 16). Since we focus on the performance 
comparison between the current PPS technique [5][8] and 
the proposed one, the initial pilot position is assumed to be 
perfectly detected at the receiver. Note that this assumption is 
reasonable due to the guarantee of the present pilot detection 
techniques in [8]-[ 10]. To make a fair comparison, the same 
values of SNR, which is denned as E^jal, have been used. 

Fig. 3 depicts the BER performance of two PPS techniques 
in an AWGN channel. At low SNRs, the BER is almost 
determined by noise and in consequence the superiority of the 
proposed PPS technique has not emerged. However, with the 
increase of SNR, the effect of the deviation of the constellation 
caused by the distortion becomes dominant; the performance 
of the proposed PPS technique is as expected overwhelming 
that of the current one. In addition, for a higher SNR where the 
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Fig. 5. BER performance comparison between the current PPS and proposed 
PPS techniques in a frequency-selective Rayleigh fading channel. 



noise effect can be ignored, the error probabilities of bits, ii 
and (72, block the decrease of BER because of the current PPS 
criterion, and therefore result in an error floor. Fortunately, this 
error floor is completely removed by the use of the proposed 
PPS technique. 

Fig. 4 shows the BER performance of two PPS techniques in 
a Rayleigh flat-fading channel, where perfect CSI is assumed 
to be available at the receiver. As can be seen, by jointly 
considering the effect of channel in selecting pilot position, the 
proposed PPS technique increases the reliability of the initial 
estimates compared with the current one, and thus improves 
the overall system performance. 

Fig. 5 illustrates the BER performance of two PPS tech- 
niques in a frequency-selective fading channel which is mod- 
eled by a tapped-delay line filter. In the simulation, CFRs are 
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estimated through the superimposed pilots (Chu sequences) 
and trigonometric interpolation. An exponential power delay 
profile is used with a root-mean-square (RMS) delay spread 
of 10 taps, where each tap follows a Rayleigh distribution. 
It can be observed from Fig. 5 that except for larger BER 
degradation compared with the flat fading channel case, the 
performance of the proposed PPS technique is still superior 
to that of the current one, which verifies the feasibility of 
the use of the lower bound on average BER in selecting the 
initial pilot position for frequency-selective fading channels. 
Actually, since the lower bound is tight enough that the 
approximated average BER can be roughly regarded to be 
equal to the exact average BER multiplied by a constant 
scale factor, the relative size among the approximated average 
BERs for different initial pilot positions would be the same as 
that among the exact average BERs. Therefore, based on the 
proposed PPS criterion which only cares about the minimal 
value of the average BERs, an accurate selection of the initial 
pilot position from the comparison of the approximated BERs 
is still guaranteed. 

It can be observed from Figs. 3-5 that for different chan- 
nel conditions, i.e., AWGN channels, flat fading channels, 
and frequency-selective fading channels, the proposed PPS 
technique with no iteration always outperforms the current 
one with one iteration at high SNRs. This implies that the 
complexity of receiver can be reduced by adopting the pro- 
posed PPS technique since one iteration additionally occupies 
0(N log N) computation time. On the other hand, by taking 
the flat fading channel case as an example, it can be readily 
found that the extra majority of computation in selecting pilot 
position at the transmitter for the proposed PPS technique with 
respect to the current one is the construction of {x-'m}m=o> 
which also requires about 0(N log N) computation time. 
Therefore, for higher SNRs the proposed PPS technique would 
express similar overall system complexity to the current one 
while still preserves the performance gain. 

V. Conclusions 

In this paper, we have verified the suboptimal property of 
the conventional PPS technique and consequently proposed an 
optimal PPS technique for a SC-FDE system with the FDPMT 
when the ZF equalization is adopted. The proposed PPS tech- 
nique is based on the minimization of the average BER whose 
closed-form expression has been derived. Simulation results 
have shown that the proposed PPS technique outperforms the 
conventional one in terms of the BER performance while 
having similar overall system complexity. By following the 
same principle, the proposed PPS technique can be readily 
applied to a scenario when the MMSE equalization is adopted. 
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